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We present a method for generating numerically a one-dimensional random surface, defined by 
the equation X3 = C(2^i)i that suppresses single-scattering processes in the scattering of light from 
it within a specified range of scattering angles. Rigorous numerical calculations of the scattering of 
light from surfaces generated by this approach show that the single-scattering contribution to the 
mean scattered intensity is indeed suppressed within that range of angles. 
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In theoretical and experimental studies of multiple-scattering effects in the scattering of light from randomly rough 
surfaces it is often desirable to be able to suppress the contribution to the mean scattered intensity from single- 
scattering processes: effects such as enhanced backscattering or the presence of satellite peaks become more readily 
observable in the absence of the background provided by single-scattering processes. 

In theoretical studies it is possible to separate the contribution from single-scattering processes to the mean intensity 
of the light scattered incoherently from the contribution from multiple-scattering processes. This separation of single- 
and multiple-scattering contributions is particularly simple in the context of small-amplitude perturbation theory^^^ 
and not much more difficult in a computer simulation approach^^'. However, it is not so easy to achieve experimentally. 
In the case of the scattering of light from two-dimensional random surfaces, the in-plane, cross-polarized scattering of 
p-polarized light suppresses the single-scattering contribution to the mean intensity of the light scattered incoherently. 
In the case of the scattering of light incident normally on a weakly rough one-dimensional random metal surface, the use 
of a surface whose roughness is characterized by a power spectrum (/(Ifcl) that vanishes identically for |fc| < kmin < '-^/c 
eliminates the contribution from single-scattering processes to the mean intensity of the incoherent component of the 
scattered light for scattering angles smaller in magnitude than sm~^ {ckmin/kuj) However, such surfaces are 
difficult to fabricate. 

In this paper we explore a different approach to the design of random surfaces that suppress the single-scattering 
contribution to the incoherent component of the light scattered from them, that is not restricted to weakly rough 
surfaces, and that appear to be easier to fabricate than surfaces characterized by a West-O'Donnell power spectrum. 

To motivate this approach, let us consider the scattering of an s-polarized plane wave of frequency to from a one- 
dimensional, perfectly conducting, random surface, when the plane of incidence is perpendicular to the generators of 
the surface. We recall that if the inhomogeneous Fredholm equation for the normal derivative of the single nonzero 
component of the electric field in the vacuum, evaluated on the surface, is solved by iteration, the inhomogeneous 
term yields the Kirchhoff approximation to the mean scattered intensity, a single-scattering approximation, the first 
iterate yields the pure double-scattering contribution, and so on*^^) . Consequently, if a surface can be designed with 
the property that the Kirchhoff approximation to the mean intensity of the light scattered from it vanishes for the 
scattering angle 9s in the interval {—9m,0m), all the scattered intensity within this range of scattering angles will 
be due to multiple-scattering processes. Consequently, our aim is to design a one-dimensional, perfectly conducting 
random surface for which the Kirchhoff approximation to the mean differential reflection coefficient vanishes identically 
for 9s in the interval {—9m, 9m)- The analysis required is simplified significantly by working in the geometrical optics 
limit of the Kirchhoff approximation. However, the results obtained still display the behavior sought. 

Thus, we consider a one-dimensional, randomly rough, perfectly conducting surface defined by the equation = 
C{xi), that is illuminated by an s-polarized plane wave of frequency w. The surface profile function C{xi) is written 
in the form 

oo 

C{xi)= J2 ces{xi-m), (1) 

where the {ce} are independent, positive, random deviates, 6 is a characteristic length, and the function Ci^i) is 
defined by'^' 

0, xi < -(to + 1)6, 

— (to -|- l)bh — hxi, —(to + 1)6 < xi < — m6, 
s{xi) ~ ^ —6, — to6 < a;i < to6, (2) 

— (to -|- l)bh + hxi, mb < xi < {m+ 1)6, 
0, (m-M)6<a;i, 

where to is a positive integer. Due to the positivity of the coefficient ce, its probability density function (pdf) 
/(t) = (<^(7 ~ c^)) is nonzero only for 7 > 0. 

It has been shown that for the random surfaces defined by Eqs. (1) and (2) the mean differential reflection coefficient 
in the geometrical optics limit of the Kirchhoff approximation is given by'*^' 



/dR\ _ 1 [1 + <M0o + 0,}]^ 
\d9s / 2hcos9o{cos9o + cos9s)^ 



^ siu 6^11 - sin gg) \ _|_ J - '^0 



/i(cos ^0 + cos 9s) J V h{cos 9q + cos 9s) 



(3) 



where ^0 and 9s are the angles of incidence and scattering, measured counterclockwise and clockwise from the normal 
to the mean scattering surface, respectively. This result shows that {dR/d9s) is given in terms of the pdf of the 
coefficient q and is independent of the wavelength of the incident light. It simplifies greatly in the case of normal 
incidence (^0 = 0°), 
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^ + tan2 .^^^ /(-itan^) + /(itanf) 



dej \ 2 J Ah 

(4) 

and we will restrict ourselves to this case in what follows. From Eq. (4) we find that if we wish {dR/dOs) to have the 
form, say, 

we must choose for 7(7) 

. / 0, < 7 < 7„, 

where 7^ = \t&n{9m/2)]/h. From this form for f{^) a long sequence of {eg} can be generated, e.g. by the rejection 
method^^\ and the surface profile function generated by the use of Eqs. (1) and (2). 

The surface profile functions Ci^'i) generated in this way arc not zero-mean Gaussian random processes, and are 
not stationary. Indeed, the mean square height of the surface, 6"^ = (C^(xi)) — is a periodic function of xi 

with a period 26, and for m = 1 is given by S"^ = [(c^) — (c)^]ft.^6^[l + {xi/b)'^] for — 6 < xi < b. The average of this 
function over a period, 6'^^ = [(c^) — (c}^]4ft,^6^/3, can be used to estimate the rms height of the surface. Similarly, 
the mean square slope of the surface is given by = ((C'(a;i))^) — (C'(a;i))^ = [(c^) — (c)^]/i^, from which the rms 
slope can be determined. The averages (c) and (c^) appearing in these expressions, the first two moments of /(7), 
are given by 




1 + cosM^ (7b) 



An example of a surface generated in this way is presented in Fig. 1. The pdf /(7) used in its generation is the 
one defined by Eq. (6), with 6rn = 40.1°. The parameters entering the definition of the function s{xi) are b = 3A, 
m = 3, and h = 0.2. For these values of the parameters we find that Sav = 0.9A, and ,s = 0.57, so that the surface is 
moderately rough. It was sampled at the points Xp = [{p + ^)b]/N, where p — 0, ±1, ±2, . . . and N — 100, and is seen 
to consist of a succession of triangular peaks and valleys. 

To show that this random surface suppresses single-scattering processes for 161,51 < 40.1°, we have plotted in Fig. 
2 the contribution to the mean differential reflection coefficient from the incoherent component of the scattered 
light, {dR/d6s)incoh, for scattering from this surface, calculated by a computer simulation approach in the Kirchhoff 
approximation^^), with and without invoking the geometrical optics limit of the latter. The results for a total of 2000 
realizations of the surface were used in carrying out the ensemble average required for obtaining the mean differential 
reflection coefficient. It is seen that in the geometrical optics limit of the Kirchhoff approximation {dR/d9s)incoh 
vanishes for \9s\ < 40.1°. In the Kirchhoff approximation {dR/d9s)incoh is not identically zero in this region of 
scattering angles, but is quite small. The difference between these two results shows how well the geometrical optics 
limit of the Kirchhoff approximation reproduces the result of the Kirchhoff approximation itself. In this figure we 
have also plotted the total contribution to the mean differential reflection coefficient from the incoherent component 
of the scattered light, including all multiple-scattering contributions. This result for {dR/d9s)incoh was calculated 
exactly by a computer simulation approach^^-* , in which the results for 2000 realizations of the surface were averaged. 
We see from this figure that there is now a low background for \9s\ < 40.1°, due to multiple-scattering, on which is 
superimposed an enhanced backscattoring peak in the retroreflcction direction (9s = 0°), whose height is nearly twice 
that of the background at its position. The latter is the result expected when the contribution from single-scattering 
processes has been subtracted . 

Although the theory underlying the approach to generating random surfaces that suppress single scattering pre- 
sented here was based on the assumption that the scattering surface is perfectly conducting, the resulting approach 
also works very well for finitely conducting surfaces. In Fig. 3 we have plotted a rigorous computer simulation re- 
sult for {dR/d9s)incoh in the case that s-polarized light of wavelength A = 612.7 nm is incident normally on a one 
dimensional random silver surface (e(a;) = —17.2 -|- «0.498). The surface is characterized by the parameters b = 3A, 
m = 3,h = 0.2, and 6m = 40.1°. Results for 2000 realizations of the surface were averaged in obtaining this figure. 
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The strong suppression of {dR/d9s)incoh in the interval \9s\ < 40.1° is clearly seen, and an enhanced backscattering 
peak at 6s = 0° rises to about twice the height of the background at its position. 

In this letter we have presented a method for generating numerically a one-dimensional random surface profile 
function C{xi) that has the property that it suppresses the single scattering of s-polarized light from it, a property 
that in the case of a perfectly conducting surface is independent of the wavelength of the incident light. The extension 
of the present approach to the generation of one-dimensional random surfaces that suppress the single scattering of 
p-polarized light is straightforward. The method described is not restricted to the generation of weakly rough surfaces. 
Surfaces defined by Eqs. (1) and (2), with a different form of the pdf fij), have been fabricated successfully in the 
laboratory^^^ . and their fabrication appears to be simpler than that of surfaces characterized by a Wcst-O'Donnell^''^ 
power spectrum. The approach to the design of two-dimensional random Dirichlet surfaces that act as band-limited 
uniform diffusers developed recently^''^ , can be used to design two-dimensional random Dirichlet surfaces that suppress 
single-scattering processes. These and other applications of the approach described here will be presented elsewhere. 
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FIG. 1. A one-dimensional random surface profile function C{xi) obtained from Eq. (1) by the use of the pdf given by Eq. 
(6) together with a function s{xi) defined by Eq. (2) with b = 3A, m = 3,h = 0.2, and Om = 40.1°. 

FIG. 2. {dR/d6s)incoh calculated by a computer simulation approachfor the random surface displayed in Fig. 1, when 

s-polarized light of wavelength A is incident normally on it. ( ) The geometrical optics limit of the Kirchhoff approximation; 

( ) the Kirchhoff approximation; ( ) the result with all multiple-scattering contributions included. 



FIG. 3. {dR/d9s)incoh calculated by a computer simulation approach for the case that s-polarized light of wavelength 
A = 612.7 nm is incident normally on a one-dimensional random silver surface. The arrows indicate the positions of the angles 
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